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Abstract 

A plasma instability that strongly influences 
the efficiency and lifetime of electromagnetic 
plasma accelerators has been quantitatively 
measured . Experimental measurements of 
dispersion relations (wave phase velocities), 
spatial growth rates and stability boundaries are 
reported. The measured critical wave parameters 
are in excellent agreement with theoretical 
instability boundary predictions . The instability 
is current driven and affects a wide spectrum of 
longitudinal (electrostatic ) oscillations . Current 
driven instabilities, which are intrinsic to the 
high-current-carrying magnetized plasma of the 
MPD accelerator, were investigated with a kinetic 
theoretical model based on first principles. 
Analytical limits of the appropriate dispersion 
relation yield unstable ion acoustic waves for 
T j/T s << *1 and electron acoustic waves for Tj/T e >> 

1. The resulting set of non-linear equations for 
the case of Tj/T e - 1, of most interest to the MPD 
thruster Plasma Wave Experiment, was 
numerically solved to yield a multi-parameter set 
of stability boundaries. Under certain conditions, 
marginally stable waves travelling almost 
perpendicular to the magnetic field would travel at 
a velocity equal to that of the electron current. 
Such waves were termed "current waves". Unstable 
current waves near the upper stability boundary 
were observed experimentally and are in 
accordance with theoretical predictions. This 
provides the first unambiguous proof of the 
existence of such instabilities in electromagnetic 
plasma accelerators. 


I. Introduction 


The electromagnetic acceleration of plasmas 
is a process that is most prone to current driven 
instabilities. 

In accelerating a plasma by electromagnetic 
means, kinetic energy is invested in the electrons 
from tne electromagnetic field. This energy is then 
imparted to the plasma by means of collisions. 
Inelastic exciting and ionizing collisions with 
neutrals are a necessary expenditure of energy that 
is needed to produce the plasma and ideally must be 
recovered as directed kinetic energy along with any 
excited internal and random translational modes of 
all heavy species through dethermalization (e.g. 
collisions with a nozzle wall or streaming along 
the lines of an applied magnetic nozzle) to achieve 
high efficiency. It a substantial recovery is made, 
we can talk of a hybrid 
electromagnetic/electrothermal accelerator. 
Dethermalization in these high power devices 
where the plasma is highly ionized is complicated 
by the associated non-equilibrium rate processes. 
The realization of a hybrid accelerator is a current 
research problem. On the other hand, for a purely 
electromagnetic device it would be tolerable to 
accept that all random and internal modes stay 
frozen, if thermalization is kept .to a minimum i.e 
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one would like to pay just the ionization price with 
the smallest possible expenditure in heating. If 
this was tne case, a high efficiency purely 
electromagnetic accelerator would indeed be 
possible since at high power operation both the 
ionization and joule heating (classical binary 
collisional heating) expenditures are modest 
compared to the useful work done by the 
electromagnetic force. In practice, however, it is 
invariably found that the efficiency at the highest 
powers does not exceed 40% indicating that a 
major dissipative process other than collisional 
heating is m operation. This discrepancy is not 
much alleviated even if electrode power losses are 
taken into account. Several related anomalies 
become apparent when experimental data obtained 
with high power electromagnetic accelerators are 
reviewed (see refs. 1, 2 and 3 for instance) 

1. considerably higher ion temperature than 
would be expected from purely collisional 
thermalization (especially near the anode) 

2. electron temperature weakly dependent on 
operating conditions or geometry (2 to 3 eV) 

3. high Ionization rates that cannot be 
accounted for by classical thermal ionization. 

4. presence of doubly ionized species with 
electron temperature not exceeding 3 ev 

5. measurements of radiation losses 
(primarily line-line) substantially higher than 
what could be accounted for on classical grounds 
(ref. 2 ) 

6. presence of thin ionization fronts across 
which strong jumps in ionization levels occur 

7. presence of suprathermal electrons (ref. 2) 


8. enhanced component heating, evaporation 
and erosion. 

9. pronounced peaks in local electric fields 
fluctuation spectra 

10. a drastic decrease in the plasma 
conductivity from classical levels 

All the above phenomena can be shown to be 
the manifestation ot plasma turbulence (nonlinear 
plasma instabilities) (refs. 3 and 4) . In these 
references the nonlinear wave theories start from 
the assumption that the instability is 
current-driven since these are known to be the 
most explosive and to lead to enhanced ionization, 
strong turbulent heating and anomalous resistivity. 
It is calculated from these theories, for instance, 
that the fraction of energy that goes into turbulent 
heating of ions (raising Tj) is a considerable 
portion of the accelerator input energy. It is also 
demonstrated that the nonlinear instabilities are in 
full effect even at low values of the total current 
J (for a fixed mass flow rate) and that the very 
rocesses of ionization and heating are dominated 
y the turbulence. When the accelerator’s plasma 
becomes effectively fully ionized, the 
disappearance of wave damping collisions with 
neutrals frees the low frequency part of the 
spectrum (typically 200 to 500 khz tor an argon 
plasma) to be current-driven unstable thus leading 
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to an accentuation in the manifestation of the 
anomalies listed above. This would happen in the 
neighborhood of \ = 1 (with 5 * J/Jfi where Jfj is the 
total current through the accelerator at the 
condition of full ionization and is definnd in refs. 5 

and 6). 

If one compares the typical current velocities 
reached by the electrons for nominal accelerator 
conditions, to the ion thermal and natural 
velocities (c.f. figure 12.C) one soon realizes that 
an' instability free electromagnetic plasma 
accelerator is highly unrealistic. This fact has 
long been recognized by soviet plasma physicists 
working on the MPD accelerator. Arefev (ref. 4) , 
for instance emphasizes that the notion of 
non-turbulent electromagnetic plasma acceleration 
is contradictory. 

Since these instabilities have a strong 
impact on both efficiency and lifetime of 
accelerators they are the subject of detailed 
studies at Princeton. 

In this paper, we demonstrate the existence 
of the current driven instability mechanism 
through data obtained from an experiment that was 
specifically staged to measure unstable waves in 
their linear stage of growth. A kinetic model is 
used to interpret the data and identify the wave 
modes and instability mechanism. 


II.Theory 


where the subscripts 1 and z denote components 
perpendicular and parallel to the magnetic field 
vector which is chosen arbitrarily to be aligned 
with the z axis . F s is normalized such that: 

(3) 

Moreover, the current density J(x,t) and the 
charge density p(x.t) are related to the distribution 
functions by the following : 

j (x,t) = S q s | f 5 (x,r,t) rd J T (4) 


=> s q s n,(x,t) = l q,J f s (x,T,t} T d 1 T (5) 

All quantities with spatial and temporal 
dependences.E, B, f«, J and p are perturbed about 
their steady-state values (superscripted with o) by 
wave-like quantities (superseded with a prime) 
and represented as Fourier-Laplace transforms 
such that for a perturbation A’ of a generic quantity 
A we have : 


A'(x,t)=|e ikx d I k|-5 2 l^_A'(k 1 u)d U ( 6 ) 

where k is the wave-vector and o> is the frequency 
Either to or k is a complex quantity : 


II. a. Mathematical Model 

In ref. 6 we treated the case of linear 
electrostatic waves in a current-carrying 
homogeneous and unmagnetized collisional plasma. 
At that time no clear experimental data on the 
nature of waves and instabilities in the MPO 
accelerator plasma were at hand. The Plasma Wave 
Experiment (PWX) described in section 111 has 
since yielded a quantitative database concerning 
waves and instabilities which helped refine the 

assumptions of the theoretical model.' The 
assumptions made in the theory discused below are 
justified for a large range of the experimentally 
measured waves. Most importantly, the experiment 
was staged such that the effects of 

non-homogeneities and collisions with heavy 
species, could indeed be neglected while the 

magnetic field,, neglected in the previous theory, is 
retained in the present one. The effects of Coulomb 
collisions are not discussed in this paper so that 
the parametric dependences of the current driven 
instabilities are not clouded. They are treated in 
ref. ( 3 ). 

The evolution of the particle distribution 
function f$(x,v,t) of species s in the 

six-dimensional phase space (x,v) is assumed to 
obey the Vlassov equation : 



+ ▼ 



4J-[e(*,o + y.b] 


« s (*.▼,*) 

dr 


0 


(1) 


k = k f + i* 

(7) 

U =5 u/ f + 17 

( 8 ) 


and y and x are the temporal and spatial growth 
rates respectively. In this nomenclature unstable 
(growing) waves correspond to a positive \ and 
negative x. The Maxwell equations, linearized for 
small perturbations, relate tne wave quantities in 
the following way: 

k * E =* w B ( 9 ) 


ik « B =s/ioj'-iu/E ( 10 ) 

k- B' =0 ( 11 ) 


ik • E' » t 0 p' ( 12 ) 


where p 0 and e 0 are the permeability and 
permittivity of free space. For the special case of 
longitudinal (electrostatic) waves, we have KxE'-O 
ana the above equations can be manipulated to yield 
(see Stix ref. 7 for instance) 


1 + i — 2L. e 

1 k J a = 


<*) 

(" - n ‘ J crV 2 ) 


where m s , q g are the mass and charge of a Darticle 
otspecies s and B(x,t) and E(x,t) are the magnetic 
and electric fields respectively. The steady-state 
and homogeneous distribution function is denoted 
by f s ° and defined as: 




(2) 


r ^ 

h s 





(13) 


where cops and u^s are the plasma and cyclotron 
frequencies of species s defined below, J n is the 
Bessel function of the first kind of order n and p s 
a non-dimensional parameter that is a measure of 
how magnetized the particles are in the wave : 
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,, -f q i }' q s B° 

Ps l 'ora, J : u « = -ir- 


We shall also allow the electrons, which 
carry the current in the MPD accelerator to have a 
(14) so*hat- dnft ve,ocity u de w ‘t h respect to the ions 




K 

"~2 


(15) 


Vt s f 2T S 1 

Vt >“ [— sf ~] 
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where we use the cyclotron radius r cs , the thermal 
velocity Vf ? and the temperature T s of species s. 
Equation (T3) relates the complex frequency w of 
electrostatic oscillations to their wavenumber k 
(or wavelength). The dispersion relation, describes 
how a packet of waves is dispersed in the plasma 
as each of the components travels with a phase 
velocity v^ - co/k. This dispersion relation reduces 
to that of an unmagnetized plasma for k . —0 (i e 
propagation along the magnetic field)! For 
propagat'on exactly perpendicular to the magnetic 
held (k 2 -0), the dispersion relation of ion or 
electron Bernstein waves is recovered. When 
solved for a real k the dispersion relation yields 
tne frequency of the oscillations as well as their 
growth or damping rate. For a two-component 
plasma we can rewrite the dispersion relation 
above as a sum of the dielectric responses %s 
representing the contributions of free space, ions 
and electrons (with subscripts i and e referring to 
ions and electrons): 


j = -* n e”de 


( 21 ) 


We shall neglect the ion drift (or flow) 
hft °fnni’ ,n K t | he * C,ispe n Sion and the formalism would 

c®.erfo P ^is ble sa«sf^d:° W,n9 P ' aSma if the ,ollowin 9 

( 22 ) 


» v. 


M. I 

Tmn.? V ! 'f J h ? f L ow ve'ocity- This criterion is 
Unri/r s .? ,,s,ied m ,h . e experiment described below 5 
Under the assumptions of isotropic. Maxwellian 

PtaSm “ *' “ d bec °™ 


Xi (u» = Jii_[i + f . 2(Ci) ] 


(23) 


k 1 u* 


X e (u,k) - _i_ _£e_ _i_ [ j + e v< lo(ih) z 


(C e ) 




1 + X e + Xj = 0 


124) 

(16) where L, is the modified Bessel function of the 
first kind of order n and of a generally complex 


. We now consider onlv .u * com P |e * parameters C, and are 

satfsfy the following criterion: frequencies th at p ^f 0 s c fj ti g| l . oc,ties n on-dimensionalized by thlrmal 


u . d << w « w ce 


(17) 


c. = 

1 ^7 


iad 


<e = 


u - k*a 


l de 


z te 


(25) 


a condition shown in section III to be satisfied by 
the experiment. In other words we consider the 

!°"f „ t0 - 1 ?® unmagnetized in the waves, which is true niaemt P 8v^ r a * the , ar 3 u ments of 2(; s ) called the 
when either: plasma dispersion function defined ft S the complex 
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(26) 


or (w r . 7) » 


(18b) 


contribution^ equations ‘(13) ® a nd° j 1 ? 6 j °b e c o m Is : i0n 


Xj(". k) = ■ 


±1 Ul 1 d»T 

u — k-r 


(19) 


The summation in equation (24) is over the 
electron cyclotron harmonics all of which except 
the zeroth order one can be nealected if : 


n de S T te 


(27) 


.. shall furthermore specialize the 

Version relation to our particular case by first 
oartfr rH & c / h 3 !- th 1 steady-state equilibrium 
&axwellian: bU On functl0ns are isotropic and 


a condition also respected by the experiment The 
dispersion equation (16) can thus be finally writter 


X + o e [1 + e^OOc, Z((.)] + Q .[i + (.z ((.)J = o 


F , K' **) = (- 


pj- ( T 1 + r z) where the parameters a s are defined 


(28) 


( 20 ) 


“s = k J 


i _ i 

< Me k \’ s k, ^ s 


(29) 


i, - 
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and are a measure of the wavelengths compared to 
the Debye lengths In our experiment both c^p 

and a,* are large on the order of 10 s and 1CK 
respectively. The major assumptions made so far 
along with the criteria for their justification are: 


1. linear waves (see equation 68) 

2. uniform plasma: 


^ ^3 d In no 

^7 q s B° 


md I VI >>•* -i jrc B .° Jj 4 In no 
11 dx ^ dx 

(31) 

3. electrostatic oscillations: -jp- << c 

(32) 

4. steady-state background: = 0 

5. plasma is collisionless in the waves: 

(33) 

** « « 

(34) 


( 

6. ions unmagnetized in the waves: 


u - k-n 


K.I-I-T7 


de 


>> 1 


i te 


(38) 


Md Ci = << 1 

* ti 


(39) 


w . '• as ,y m Ptotic expansions of the olasma 
dispersion function Z(C~) t or lame sn /S 
arguments can therefore be used : 96 3 d S a 



If, I » 1 


(40) 


z (fj) * 2C S + -j- (j + i/T " exp fj, |f s | << I 


(41) 

Moreover we study only the case of 
long-wavelengths or equivalently the case of 
electrons strongly magnetized in the wave i.e. u e << 
1. For small arguments the modified Bessel 
function of the first kind and of the zeroth order 
can be approximated bv: 


<< k v, 


ti or u d <<: ( u c 7) 


( 35 ) 


«^I,(P,) = 1-fle for It, « l 


(42) 


7. electron cyclotron harmonics are smeared: 


“de 5 * te 

than the (Tow® PhaS %£dg? es consi(jerab| y lar 9 er 

-Z- » r, < 37 


The dispersion equation can now be solved to 
arrive at a parametric description of waves phase 
velocities, criteria for their instability and their 
rates - . Analytical solutions of equation 
(28) are possible in the limits of very large and 
very small values of Tj/Tp These two limiting 
cases are reviewed first. THe first case leads to 
electron acoustic waves and the second gives ion 
acoustic waves. Both of these modes can oe made 
unstable by the electron current. The magnitude of 
Tj/T e that is necessary for the existence of each of 
these two modes is expressed in terms of the 

plasma and wave parameters. The plasma state and 
probe geometry of the experiment reported in this 
paper were such that neither of these two modes 
could be observed. Therefore the existence of 

current-driven instabilities could not be 
ascertained using these analytical solutions. This 
necessitated the treatment of the case of moderate 
Tj/Tp. Numerical solutions of equation (28) in the 
complex domain were thus undertaken. In this 
paper, for the sake of clarity, we shall only discuss 
the solutions near and at the instability boundaries 
which fortunately are sufficient to interpret the 

experimental data and demonstrate the existence 
of current-driven instabilities. For the complete 
solutions of equation (28) including substantial 
growth and comparison to experiments we refer the 
reader to Choueiri (ref. 3 ). 


U-& Unstable Electron Acoustic Waves 
(T|/T e >>1) 

We first assume that T\/T e is sufficiently 
large so that : 7 


eiecirons, equations (23 
respectively : ; 
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z ce 
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k 3 Z> 2 
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where 

u = w-k-u^ 


(43) 


(44) 


(45) 




‘ v ie A cos ® + ‘TOrfg co 


( 46 ) 


7 = -v j e T ^ 1 kcosfl | [ j-. p l. j e~ A*Ti/T, 


T ti 
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bi&w S a innn e ^K°.K acous,ic velocity and is defined 
oeiow along with the parameter A 




- 1/2 
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0 is the angle between the wave vector k and the 
magnetic field while <t> is the angle between k and 
u de, (see figure 11 for the case of the M PD 
accelerator), o is the sign of the dot product k.ud- 
• The instability boundary can be cast in terms ofa 
threshold current velocity by setting v to zero in 
the above equation to obtain: ' ' 

•L-h (-^-)' , ’'- T,A ’ /T --«-']- i V 

(49) 

The above equations indicate that electron 
acoustic waves cannot be driven unstable (i.e. r>0) 
by the current unless o<0. Although these waves 

m flfr r A S h }5 the i^ M u D P ,asma ]since their 
instability threshold could be satisfied especially 
near the anode where the ions are hot ) they are not 
relevant to the experimental case reported here 
since wave injection in PWX was always in the 
same sense as the current i.e. a was always 

Prnm A 0 ;af a Eq / at i On o S < 46 u and 47 ) are recoverable 
from Arefev (ref. 8 ) when assumptions (31) are 
made. Finally we recast the criterion T;/T e >>1 in 
more specific terms using (38 and 46)* 


Finally, we recast the conditions (50) in 
closed form using equation (52) to give us more 
explicit criteria for the existence of these modes 


cos0 > 


(4H 


1/2 


<< l + kV 


( 56 ) 


(57) 


T. 


>> 1 + 


k 3 cos 3 0 


(50) 


II. b Unstable ion acoustic waves (Tj/T e <<1) 

T ih, parallel to the above case we assume that 
J l, e. IS sufficiently small so that we have the 
following conditions 

» 1 l( e l « 1 (51) 

,,,, Using the asymptotic expansions (40) and 
(41) and proceeding as in the case above we obtain 
after some manipulations the phase velocity and 
growth rate for the more familiar Ion acoustic 


It is clear that in both limitinn 

c2r U r en? n vJn 4? r) ?" d (5 ^ > the threshold on the 
acr^ofatn 0 r 0Cl h y ,s easily surpassed in the MPD 
? he L e cur ( ent velocities are. in many 
inn ^h 5, tv li 0 °, rd ? rs of rnagnitude larger than the 

susceotibmtv Ve of CI, th B A ni hQugh ‘ he y emphasize the 
! n .i.kr.“ y ,he P'asma to current driven 

equaMon 0 f S 28l ,h r a nnn? V h analy j' cal solutions of 
Dw/v Ca - nno used to interpret the 
fill PW ? s,n 5 e .. th ® inequalities in (50) and 
re not _, satisfied for the experimental 
conditions under which the instabilities were 

parameters) 0 vfc"!J 2 rif , for 5 lasma and P'°be 
fh« a ir eters ^ e are therefore driven to consider 
the more complicated case of moderate Tj/Tg 


II. c The case of moderate Tj/T e and the 
definition of "current waves" 

We first rewrite equation (28) in more 
convement form by noting tnat for most situations 
2L2 ler «* t0 the experimentalist the measured 
wavelengths are considerably larger than the Debye 
iSE? « > From the - deflnition at . (29) this translates 
I- 1 ?®. earl ’ er * this condition is 

S^£ Sa # S ^A by yhere a e and ctj are on the 
orders of 10° and 10' respectively. Dividing 
equation (28) by aj with aj » Iwe obtain from the 
real and imaginary parts the following set of 
nonlinear equations relating the complex variables 
Ce and £i with Tj/T e and ji 9 as parameters: 


y si A 


( 52 ) 


T 

1 + [ 1 + e^Io0*K e Re { z« e ) 

+ <i Re ( Z(Cj) ) - 0 


7. to, 

e 1 ,J jz(g) + C i Im(z(C i )]=0 


( 58 ) 


+ ] a - u de cosS . * J 1 


( 53 ) 


btlowf along i °w r ith a A'° US,iC Ve '° city and is de,ined 


r t.i 1 / 2 

0 r t 

[-:) • 

A -[l + rt*,] 


- 1/2 


( 54 ) 


Even though the growth or damping rate 
depends on the sense of propagation a (as seen from 
equation 53), the threshold for instability is 
independent of o, unlike the case of electron 
acoustic waves. This is readily seen by setting y - 
0 and solving for the critical current velocity at 
the boundary ’ 


i 1 / 2 r T - I 3 / 2 „ 


■i-’.["*'-Rr HH 


A T e /Tt 


( 55 ) 


In order to numerically solve the above 
equations we need the following numerical tools. 
First, a plasma dispersion function generator in the 
complex domain for which the continued fraction 
representation of integral (26) was used (ref. 9 ). 
For J lm(C s ) | < 1 this method is weakly convergent 
and the Lanczos-x method of Franklin (ref. 10) was 
used instead. Second, the modified Bessel function 
of the first kind and of order zero with a generally 
complex argument was generated by means of an 
ascending series for small to moderate arguments 
and an asymptotic expansion for large arguments 
Finally in order to find the roots of the nonlinear 
set of equations a specialized recursive 
Newton-Raphson algorithm was developed. 
Reference 3 has more details on these numerical 
tools. 

The extensive solutions of these equations 
for substantial growth or damping can be found in 
ref. 3 . In this paper, we present the more 

manageable set of solutions at the instability 
boundaries which is sufficient to identify the 
measured instability. For a marginally stable wave 
the plasma dispersion function (26) with a real 
argument x s becomes 


Z (*,)=- e 1 [ 2 J fl e t5 dt-i(r)V 2 ] 


( 59 ) 
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. An interesting singularity happens when the 

and upon substitution in (58) we obtain the curves for xj and Ixq 1 meet. This can be seen in 

nonlinear set of equations at the instability !'9 ure ? ( 4 ) and (5) for the cases of upper and lower 

boundary ' boundaries respectively. 


l + - T i-[l-2e^ I 0 ( / i,) Xe | a 

r x i 

+ 2x. c l Xl dt = 0 
1 J o 

e^IoW x e t~ xl + * x? = 0 


whers all the variables and parameters are 
non-dimensional. The real variables are 




u - k-u 


de 


Vti 


( 60 ) 


( 61 ) 


In figure (1) the roots representing the 
non-dimensional phase velocity at the upper and 
lower limits of the instability are plotted versus 
the temperature ratio Tj/T e with n e as parameter. 



n^T* 


Figure 3. xi and |x e | at the instability boundary as 
a function of the temperature ratio for the cases of 
nil and moderately small n e . 


WSTAMITT BOUKDaAIQ 



Figure 1. xj at the boundary of instability as a 
parameter^ 100 tem P erature ra «o with p 9 as 


Figure 4. Solutions at the upper boundary showino 
current wave nodes (xi-|x e j). 


For a fixed ^i e there exists a cutoff 
temperature ratio for which the upper and lower 
boundaries merge. This singularity recedes to 
higher values or Tj/T e as the electrons become 
more magnetized in the span of a wavelength i.e. as 
diminishes. For the extreme case of p 0 - 0 the 
plot shows the lower instability boundary (for 
Tj/T e up to 6) while the corresponding upper 
boundary is at infinity. For a fixed ll 0 and for 
temperature ratios below the cutoff value, the 
curves indicate the highest and lowest phase 
velocities of the unstable oscillations. The 
existence and physical meaning of other roots are 
discussed In ref. 3 . Figure (2) is a similar plot for 
|x e |. The curve for 4 q *0 is not shown here for the 
sake of clarity instead it is shown in figure (3) 
along with the curves for a moderately high ji e . 


INSTAStUTT LQVtJt tOUMDAJUtS 
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2 -J x <?[ at * h8 boundary of instability as a 
parameter 0 tem P erature with J. as 


Figure 5. Solutions at the lower boundary showino 
current wave nodes (xj-|x e j). 


x i * I x 9 I (61) we get the phase 
velocity of the wave at the threshold of instability 


(-T-) =' 


u de 


cos^ 


• cos# 


rtr-i 


TF 


+ i 


(62) 


and for waves propagating almost perpendicular to 
the magnetic field or more explicitly at an angle 8 


costf << 



“. J 


(63) 
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their stability boundary is simply related to the 
current velocity 

(-£-)* = U de COS0 ( 64 ) 

Definition: We shall call "current waves" 
those waves that become marginally unstable when 
their phase velocity in the direction 
6<<cos' 1 (T[m e /T e m j) 1; 2 becomes equal to the 
current velocity component along the wave vector. 

Conditions for the existence of current waves 
can be found by setting xj - |x e [ in equations (60) 
(see ref. 3 for a more detailed discussion). 
Graphically, the current wave branches can be seen 
by combining the two plots of figures (4) and (5) 
and connecting the nodes as shown in figure (6). 


CURRENT WAVE BRANCHES 



Figure 6. Branches of the "current wave" mode. 


Aside from at the nodes there exists a 
branch of current waves at asymptotes where xj -► 

I x e | . In anticipation of the data of PWX where the 
experiment was staged in a nearly isothermal 
region of the plasma, the case of an isothermal 
plasma (Tj/T e - 1) was solved for and the upper 
and lower instability boundaries are shown in 
figure (7). From this plot we conclude that current 
waves are to be expected in an isothermal plasma 
for ii e 5 .1 as seen from the upper boundary curve. 
At sufficiently low ii 0 current waves are to be 
expected even if tne plasma departs from the 
isothermal condition. In PWX, p Q was typically on 
the order of 10’ 4, and indeed a current wave was 
observed and its instability measured as discussed 
in the next section. 


1N9TABIUTT BOUNDARIES TOR AN CSOTHERMAL PLASMA 



Figure 7. Boundary solutions for an isothermal 
plasma as a function of j i e showing an asymptotic 
branch of current waves at small values of n 0 . 


III. Experiment 

The goal of the experiment is to supply 
measurements of plasma dispersion relations and 
spatial growth rates under conditions compatible 
with the theory so that the wave modes and 
instability mechanisms could be identified. Two 
types of measurements were made for each 1 msec 
pulse of the accelerator and are referred to as PWX 

and PCX. In PWX (Plasma Wave Experiment), wave 
parameters were measured such as dispersion 
relations (phase velocity spectra), intensity 
spectra, spatial growth rates, etc... PCX (Plasma 
Characterization Experiment) was concerned with 
the plasma parameters such as electron 
temperature, charged species density and magnetic 
field. The execution of these measurements was 
tailored such that all the above listed parameters 
were obtained during the quasi-steady part of the 1 
msec pulse. Improvements on the hardware 
software and methods reported in ref. { are 
described below. 


III. a Experimental set-up 

The plasma accelerator selected for these 
studies was the Full-Scale Benchmark Thruster 
(FSBT) shown in figure (8). All experiments were 
conducted in a fiberglass tank 1.8 m in diameter 
and 4.8 m long maintained, by diffusion pumps at a 
pre-run pressure of a^bout 10- 5 torr. A 3000 p F 
pulse-forming network produces 1 msec long 
rectangular pulses driving a current through the 
accelerator of 3 to 30 kA. Argon is injected from 
the accelerator backplate at a mass flow rate of 6 
g/sec. The Drobes for PWX and PCX are shown in 
Figure (9). The PWX probes consisted of an emitter 
and receiver probe, each having two parallel 
cylindrical tungsten wires 5 mm long and .127 mm 
in radius. Low attenuation RG-217 coaxial cables 
were used to connect the probes to the diagnostics 
all of which were housed inside a Faraday cage. 
The current fluctuations of the probes were 
measured by Tektronix current probes and recorded 
on the two channels of a Nicolet digital recorder. 
Sampling rates of 200 nsec were selected for the 
results reported here corresponding to a Nyquist 
frequency of 5 MHz. The recorder is capable of 
storing 4096 samples. To prevent aliasing, a 3 MHz 
analog low pass filter was connected at each input. 
An analog spectrum analyzer (HP 8553b) with a 1.2 
GHz ceiling RF section was kept on line with the 
emitter probe to provide real-time monitoring of 
the frequency spectrum. The frequency response of 
ail the analog instrumentation was measured with 
the spectrum analyzer and found to be flat below 
the foldover frequency. 

Since all plasma oscillations are in a steady 
saturation state, the dispersion relation of the 
plasma can only be obtained by monitoring the 
phase shift and amplitude variations of an injected 
wave packet. This injection was done by 
connecting the emitter (upstream) probe to a 



Figure 8. The Full Scale Benchmark Thruster 
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Plasma Wave Driver which is an RF oscillator that 
can provide a harmonically rich spectrum of 
oscillations of about 20 volts peak-to-peak 
centered at the plasma potential. The receiver 

probe was held at a constant bias to draw ion 
saturation current. The PCX probe cluster 
consisted of an asymmetrical double Langmuir 
probe positioned parallel to the flow and a 
magnetic probe whose coil axis is positioned 
parallel to the local magnetic field. The Langmuir 
probe was chosen to be asymmetrical in order to 
increase the fraction of the electron distribution 
function that is sampled. The two cylindrical 
components were 7.62 mm long and .127 and .0635 
mm in diameter. The bias voltage on the double 
probe was swept from -30 to +30 volts in about 50 
psec. The sweeping is triggered 100 jisec into the 1 
msec pulse. The swept response of the current in 
the probe circuit was measured by a Tektronix 
P6041 Hall type current probe that has a frequency 
response of dc to 50 MHz. The probe circuit voltage 
was measured by a P6015 Tektronix voltage probe 
that also assured a capacitive decoupling of the 
floating Langmuir probe from the instruments. The 
measured probe characteristics were recorded on 
the two channels of another Nicolet digital 
recorder. The accelerator’s terminal 
characteristics (voltage and total current), the 
output of the magnetic probe integrator and the 
various trigger pulses were recorded on a 
multi-channel Siomation 1015 digital recorder. In 
all up to 10 digital channels were in operation for 
each pulse. Tne data were transfered to an IBM 
9000 micro computer for storage then up-loaded to 
the Princeton University IBM 3081 mainframe 
computer where all of the processing software is 
resident. Data processing ana graphical 
manipulations are monitored on an interactive IBM 
3277 GA workstation. 

The PCX software has been described in (ref. 
6) and various progress reports (ref. 11). The 
software processes the probe characteristics for 
an estimate of the electron temperature using 
asymmetrical double probe theory (ref. 12). The 
charged species density is obtained from an 
interactive computerized version of Lam's Unified 
Theory (ref. 13 ). These data, along with other 

experimental inputs, are used to calculate various 
state, collisional and natural parameters of the 
plasma at the probe’s location from established 
classical models. The probe theories validity 
parameters are then evaluated for consistency. The 
output is a list of the parameters and a set of scale 
ordering diagrams as later shown in figure (12). 



Figure 9. Probe duster for PWX and PCX. 
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The PWX software uses a Fast Fourier 
Transform algorithm to compute the phase and 
intensity spectra of the recorded oscillations. 
From these spectra and the probe cluster geometry, 
various spectra of parameters relevant to the 
studv of instabilities can be estimated (see ref. 
(19 ) an ^ 3re ^ ottec * as sbown ' n figures (13) to 

A background clipping is done for all the data 
in order to eliminate the scatter due to the 
incoherent oscillations of the background plasma. 
This is done by clipping the background in the 
frequency domain below a certain attenuation. This 
process is completely unbiased as no 
pre-determined selection of peaks is made. 
Clipping could not help reduce tne scatter when the 

g eaks of the injected waves were swamped in the 
ackground noise . This type of severe 
deterioration to the signal-to-noise ratio happens 
when the accelerator’s $ is above .75 for the case 
where the emitter probe is 2 cm downstream of the 
tip of the cathode. A more powerful wave driver is 
needed to overcome this limitation. All PWX 
instrumentation tested negative for spurious phase 
shifts, attenuation or amplification for all 
frequencies up to 10 MHz. 


IH.b Experimental Results 

results reported here were obtained at £ 
- .67 which, for the full scale Benchmark thruster 
device, corresponds to operation at 11.8 kA with 6 
a M 9 n°* n * J he P rob ® cluster position is shown 
IU*I° ur ? i 1 ^ and was such that the emitter probe 
” a L a f 1-7 cm downstream of the tip of the 
joining the two probes which 
represents the direction of the measured wave 
vector was parallel to the accelerator's axis of 
symmetry and at a small radial distance from the 
axis to avoid any r-0 singularities. The probe 
^®P?r a V on a compromise to various 

conflicting effects among them: diagnostics 

sampling capability, frequency response, bandwidth 
of unstable waves and the assumption of uniform 
plasma. The probe cluster radial location was 
chosen near the cathode axis for two main reasons: 
i) isothermal plasma conditions are most satisfied 
there as indicated by recent spectroscopic 
measurements of Tj and T e (ref. 14) and 2) the 
higher electron Hall parameters (Q e - co ce /v ce \ 
( c J’ 9 u re 10 ) permits simplification of 
on !S used to verify the wave mode 

at the boundary of instability. A high hall 
parameter implies that u y » u x in figure (11) and 
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the three vectors u d0 , k and B become coplanar. 
This combined with the fact that for a self induced 
field, B and u de are always perpendicular imply 
that (cf. figure 11) 


* * 



(65) 


and the phase velocity of current waves at the 
stability boundary becomes simply 




u de 



( 66 ) 



Figure 11. 
accelerator 


The k, u and B 
plasma. 


vectors in the 


Figure (10) is a map of the current contour in 
the plasma that was constructed from magnetic 
probe measurements (ref. 15) with the FSBT at 15.3 
kA and 6 g/sec (i.e. 5 - -97). This map can be used 
to estimate the average current velocity u d0 in the 
experimental control volume bounded axially by the 
cross sections a-a and b-b of figure (10) where 
the probes are located. Assuming that the current 
density in the control volume scales linearly with 
the total current (this is a good assumption when 
the shape of the control volume bounded by the 
probes axial positions, the line of constant 
enclosed current and the symmetry axis stays 
similar with the total current which is the case for 
4 £ 1} the current velocity at an arbitrary £ can be 
obtained from the following formula: 


u 


de ~ 


4 J o r 

5 en t TX I 
0 



(67) 


where 5 0 and J 0 are respectively the value of E and 
the magnitude of the nearest measured enclosed 
current line for the conditions under which the map 
was obtained, r a _ a and rw-b are radii to this 
line from the axis at the locations of the emitter 
and receiver probes respectively, e is the 
fundamental charge and nj is the ion number 


density measured Oy PCX for the case in question 
in this experiment we have $.= . 87 from the 
conditions ot the map. r-,- .6 SSt. ru* -12 cm 
m 1 kA from figure f 10) , £ = .Sr and n| - 
2 5 • 2 1 x 1 0 2o m- 3 from the output of PCX (figure 12 a ) 
thus yielding u de - 120.2 km/sec. h 

The experiment conditions and the plasma 
parameters are shown in figure (12) . In computing 
the various plasma parameters from the measured 
quantities, an isothermal plasma was assumed. An 
ionization fraction of 0 5 was also assumed for an 
c U the . fre 9 uencies of Collisions with 
Saf ] a value ls - 9 ") 3,1 of which where 
substantially lower than the measured wave 
frequencies. The local Dlasma flow velocity was 

estimated with three independent methods (ref. )j ) 
at 5.83 km/sec. 


Figures (13) to (18) show the results of the 
Plasma Wave Experiment. The time domain traces 
of the emitted and received signals are shown 
along with their corresponding power spectra in 
the frequency domain in figures (13) and (14) . The 
cross spectrum and the phase spectrum are shown 
in figures ( 15 and 16). The data scatter in the 
phase spectrum was eliminated in figure (16 ) by 
background clipping at *24 db indicating that the 
scatter was indeed due to the incoherent plasma 
noise. The time delay spectrum was also plotted 
and with known probe separation was used to 
arrive at the phase velocity spectrum which in turn 
is equivalent to the dispersion relation plot 
(figures 17 and 18. a respectively). Finally, the 
clipped power spectra at the emitter and receiver 
were used with the known probe separation to 
arrive at the sought spatial growth rate spectrum 
of figure (18.b). In this plot negative values of the 
ordinate (the imaginary wave number) correspond 
to unstable waves (c.f. equation 5). Waves grow or 
get damped exponentially in space as e** x where k 
is the ordinate of the spatial growth rate plots. 

The spatial growth spectrum shows that the 
local plasma is unstable for oscillations between 
1.7 and 2.7 MHz. This, to our knowledge, is the 
first direct experimental measurement of a 
convective plasma instability in electromagnetic 
plasma accelerators. 

That the instability is an inherent trait of the 
plasma and not an effect of the active wave 
injection (i.e. that the injected waves are indeed 
linear) can be easily checked using the above 
plasma and wave parameters in the following 
criterion (ref. 16) w ! 

ky ts * (Nk’A^) 2 N> ^ 

where N is the number of plasma periods up to 
which linear waves can be assumed. This in 
essence a criterion for the smallness of the 
injected wave energy density compared to the 
plasma thermal energy density. The above criterion 
is amply satisfied in this experiment. 


The strong damping at lower frequencies (cf. 
figure 18.b) is mainly due to collisional damping. 
Below 700 kHz no injected peaks survived in the 
background as can be seen from figures (14 and 15) 
All spectra were therefore clipped at that 
frequency to avoid cluttering the plots with the 
low frequency incoherence. When the plasma 
becomes fully ionized ( 5-1 ) and the damping 
effects of collisions with neutrals are eliminated, 
it is to be expected that the lower part of the 
spectrum also becomes unstable. 

Waves above the upper instability boundary 
were too damped to be observable in the narrow 
band between the boundary and the cutoff of the 3 
MHz low pass filter. The location of the boundary 
(i.e. the frequency there) was verified to be 
insensitive to the nearness of the low pass cutoff 
and was thus a true plasma effect. 
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That we should expect a current wave mode is 
dear from the above theory. As can be seen from 
figure (12.b) the condition Q) C j<<£d<<o) C0 is amply 
satisfied for our experiment (cu between 1 and 3 
MHz) and therefore the above theory can be invoked. 
The value of the parameter p e easily computed 
from the above measured quantities, is less than 
10* 3 for all measured waves (the maximum 
wavenumber being 140 m- 1 ). This means that the 
electrons are strongly magnetized on the scale of a 
wavelength and from the theory above (cf figure 7 
for instance) , the asymptotic branch of current 
waves is to be expected. The wave phase velocity 
at the instability boundary would, according to 
(66), be very close to the electron current velocity 
Ucfe which for our case is 120.2 km/sec. This is 
indeed the case as shown in figure (19) where the 
theoretical prediction of the critical phase 
velocity is indicated by a threshold line on the 
phase velocity plot and equivalently a slope on the 
dispersion relation plot. It is seen that at the 
. boundary of instability defined by the spatial 
growth data (i.e. 2.7 MHz) the phase velocity 
coincides with the predicted one. this proves that 
the measured instability is current driven. 


IV. Conclusions 

The major wave mode and instability 
mechanism inherent to the electromagnetic 
acceleration of a plasma are Identified 
experimentally and theoretically. Current 
electrons which would ideally impart momentum to 
the plasma ions by the mildly dissipative process 
of Coulomb collisions create, in addition , an 
unstable distribution of energy that fuels 
disturbances into unstable exponential growth. 
Such current driven instabilities are notorious for 
dissipating energy in strongly enhanced ionization, 
turbulent heating of the ions to excessive 
temperatures, anomalous conductivity, component 
erosion and strong oscillations in the accelerator 
plasma limiting both efficiency and lifetime more 

than any other type of instability. In this paper, a 
theoretical study based on kinetic theory was used 
to interpret experimental measurements of an 
instability in the MPD thruster plasma. 
Phenomenological traits of the excited waves near 
their stability boundaries led to the unambiguous 
identification of the mechanism behind their 
instability. 
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Figure 12. Output of PCX for the present 
experiment: a) List of plasma parameters, b) 
frequency scale ordering diagram and c) velocity 
scale ordering diagram. 
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Figure 14. Frequency domain 
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Figure 15. Unclipped power 
cross-spectra (Doppler shifted) 


Figure 18. Dispersion relation 
phase (spatial growth spectra) sh 
(Doppler shifted). 
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